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ABSTRACT: Severe acute respiratory syndrome (SARS) was a worldwide epidemic caused by a coronavirus 
that has a cysteine protease (3CL?°) essential to its life cycle. Steady-state and pre-steady-state kinetic 
methods were used with highly active 3CL?° to characterize the reaction mechanism. We show that 3CL° 
has mechanistic features common and disparate to the archetypical proteases papain and chymotrypsin. 
The kinetic mechanism for 3CL?"°-mediated ester hydrolysis, including the individual rate constants, is 
consistent with a simple double displacement mechanism. The pre-steady-state burst rate was independent 
of ester substrate concentration indicating a high commitment to catalysis. When homologous peptidic 
amide and ester substrates were compared, a series of interesting observations emerged. Despite a 2000- 
fold difference in nonenzymatic reactivity, highly related amide and ester substrates were found to have 
similar kinetic parameters in both the steady-state and pre-steady-state. Steady-state solvent isotope effect 
(SIE) studies showed an inverse SJE for the amide but not ester substrates. Evaluation of the SIE in the 
pre-steady-state revealed normal SIEs for both amide and ester burst rates. Proton inventory (PI) studies 
on amide peptide hydrolysis were consistent with two proton-transfer reactions in the transition state 
while the ester data was consistent with a single proton-transfer reaction. Finally, the pH-inactivation 
profile of 3CL’" with iodoacetamide is indicative of an ion-pair mechanism. Taken together, the data are 
consistent with a 3CL? mechanism that utilizes an “electrostatic” trigger to initiate the acylation reaction, 
a cysteine—histidine catalytic dyad ion pair, an enzyme-facilitated release of P}, and a general base- 


catalyzed deacylation reaction. 


Severe acute respiratory syndrome (SARS)! was a world- 
wide epidemic that appeared in November of 2002. The 
overall mortality rate was estimated by the World Health 
Organization to have reached 15% of all patients and up to 
50% for patients over 65 years of age. In the wake of 
significant media attention, the research community applied 
intense and productive attention to this new viral disease. 
Four months after the first SARS case was identified, a new 
coronavirus was shown to be the causative agent: SARS- 
CoV (/). Within | month, the entire genome of SARS-CoV 
had been sequenced (2, 3). The SARS-CoV was shown to 
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be phylogenetically distinct from the three previously identi- 
fied coronavirus groups (4). The main protease responsible 
for maturation of the polyproteins and central to the life- 
cycle was identified as a cysteine protease similar to the 3C 
family of cysteine proteases: 3CL?'° (also known as MP"). 
The 3CL?” has been cloned, expressed in Escherichia coli, 
and purified (5, 6). Its enzymatic activity has been measured 
by the development of both discontinuous and continuous 
enzymatic assays (6—/1). Initial kinetic studies showed that 
the 3CL?” had modest activity (Kea/Km of 177 M7! s~!) (6), 
while other reports had high protease activity (Kea/Km > 10000 
M7! s“!) (8, 12). Differences in enzymatic activity and 
dimerization constants may be attributed to the 3CL?® 
construct design, the substrate, or the handling/storage 
conditions (//, /3). 


The crystal structure of the protease has been released 
(PDB entry 1q2w by Bonnano et al. 2003) with other 
structures published subsequently (/4, 15). 3CL*° has been 
characterized as having a chymotrypsin-fold (/4, 16). The 
structure revealed that the full-length protease contains a 
catalytic domain (aal—184), a linker region (aal85—200), 
and a putative regulatory domain (aa201—303). Structural 
studies have provided some insight into the 3CL? reaction 
mechanism. The active protease is a homodimer with the 
Cys145—His41 catalytic dyad located in the cleft between 
the first two domains (7, /4—/6). Detailed studies of the 
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functions of the individual domains have revealed that the 
carboxy-terminal domain (aa201—303) has a critical role in 
the formation of the active dimer (/7). 


While at least eight superfamilies of cysteine proteases 
exist, our understanding of cysteine protease mechanism is 
derived primarily from studies of the papain family of 
cysteine proteases or inferred from related serine proteases 
(18). Papain is the prototypical cysteine protease, which 
catalyzes peptide hydrolysis with a catalytic dyad ion pair 
(Cys25—His159), which is oriented by a third residue 
(Asn175) (/8). For the acylation reaction of papain, the ion- 
pair dyad is the catalytic device responsible for the nucleo- 
philic attack on the substrate’s scissile amide carbonyl, and 
it also is involved in the facilitation of the leaving group 
stabilization (JS—2/). For papain, His159 acts as a general 
base to facilitate the activation of an active site water in the 
deacylation half-reaction (/9). With 3CL°", an analogous 
catalytic dyad (Cys145—His41) has been observed in the 
crystal structure. But unlike papain, the histidine residue of 
the 3CL?® catalytic dyad is oriented by a water molecule 
not an asparagine residue (/6). The other common enzymatic 
mechanism for proteolytic peptide cleavage is general base 
catalysis. The fundamental feature of a general base mech- 
anism is that the catalytic group participates in a proton- 
transfer that stabilizes the transition state. General base 
catalysis is the commonly accepted mechanism for chymot- 
rypsin (22). Chymotrypsin has a catalytic triad (Ser195/ 
His57/Asp102) that acts as a charge relay system to increase 
the nucleophilicity of Ser195 through the abstraction of the 
hydroxyl proton in a concerted attack of the substrate 
carbonyl of the scissile bond. 


The minimal model for cysteine protease catalysis (Scheme 
1) is a three-step reaction that first involves the binding of 
the substrate to form a Michaelis complex (ES). Next, an 
acylation reaction occurs with the simultaneous formation 
of the covalent enzyme acyl intermediate (ES’) and the 
alcohol or amine product (P}). Finally, a deacylation reaction 
occurs through the hydrolysis of the ES’ acyl intermediate 
to generate the second product, P>. In this model, the 
acylation rate is described by kz and the deacylation rate by 
k3. An alternative model (five step) has been proposed that 
has an additional conformational change to facilitate the 
release of P, (Scheme 2) (23). The (ES’*P\)* complex 
represents the enzyme conformation from which the alcohol 
or amine P; rapidly diffuses to leave (ES’). 


The enzymatic mechanism for 3CL”" has been proposed 
to utilize chymotrypsin-like general base catalysis and not a 
papain-like thiolate—imidazolium ion-pair mechanism (/0). 
To support this assignment, the authors presented multiple 
lines of evidence. They showed that the mutation of the 
catalytic cysteine to serine (C145S) was active with a 40- 
fold reduction of activity. An inverse solvent isotope effect 
(SIE) was observed for the wild-type 3CL’, while a normal 
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SIE was observed for the C145S mutant. Last, the measured 
pX, values were reported to be consistent with a general base 
mechanism. This first reported kinetic study used a 3CLP'° 
construct that had an additional N-terminal residue, which 
significantly reduced enzymatic activity. 


In this paper, we have further studied the mechanism of 
the 3CL?° protease using steady-state and pre-steady-state 
methods on amide and ester substrate hydrolysis reactions. 
The rate constants that define the kinetic mechanism of ester 
hydrolysis were characterized. Included in this report are 
solvent isotope effect and proton inventory studies utilizing 
amide and ester substrates. Although proton inventory studies 
have been widely used for serine proteases, they have not 
been widely used with cysteine proteases (20). The proton 
inventory studies coupled with pre-steady-state isotope effect 
evaluation provide a deeper level of understanding of the 
recently discovered cysteine protease 3CL?"°. We also used 
the mechanistic approach for characterizing the catalytic 
mechanism of a cysteine protease: pH dependence of 
iodoacetamide inactivation. The inactivation profile is more 
consistent with an ion-pair mechanism. Taken as a whole, 
our data are more consistent with the 3CL?"° catalysis of the 
acylation reaction occurring through an ion-pair mechanism 
and the catalysis of the deacylation reaction by a general 
base mechanism. Evidence is provided for an “electrostatic 
switch” to trigger the acylation reaction and an enzyme- 
facilitated release of the first product. 


MATERIALS AND METHODS 


Materials. Peptide 1, peptide 2, and peptide 3 were 
purchased from SynPep (Dublin, CA) and had the following 
purities: >97%, >93%, and >91%, respectively. Peptide 
4, peptide 5, and peptide 6 were purchased from CPC 
Scientific (San Jose, CA) and had the following purities: 
> 95%, >93%, and >90.1%, respectively. The following 
reagents were purchased from Sigma-Aldrich Chemical Co. 
(St Louis, Mo): HEPES, BTP bis-tris propane, MES, DTT, 
DMSO, EDTA, D;O, iodoacetamide, pNA, and pNP. Tween- 
20 was purchased from Calbiochem (San Diego, CA). 


Expression of 3CL?" The design of the 3CL? was based 
on the genome sequence of the Urbani strain SARS-CoV 
coronavirus (GenBank accession number AY278741, nucle- 
otides 9985—10902), and the synthesis of the 3CL?"? DNA 
was prepared by Blue Heron Biotechnology (Bothell, WA). 
With the resulting DNA as a template, the gene encoding 
the mature form of 3CL?® (residues 1—306) was amplified 
by PCR with forward primer 5’-GTGCTCATGAGTG- 
GTTTTAGGAAAATGGCATTCCCGTCA-3’ and reverse 
primer 5’-AAGCTTATTGGAAGGTAACACCAGAGCAT- 
TGTCTAAC-3’. The PCR products were cloned into pCR- 
Blunt If TOPO vector (Invitrogen, Carlsbad, CA). The 3CL?" 
DNA was excised from the TOPO vector and subcloned into 
the NcoI and HindIII sites of the pET28a vector (Novagen, 
Madison, WI). The resulting expression construct encoded 
the mature form of the 3CL”° (residues 1—306) with the 
addition of the initiation codon. The correct sequences were 
confirmed by DNA sequencing of both strands. 3CL?"? with 
the deletion of C-terminal two residues (residues 1—304) 
was made by substituting reverse primer 5’°-AAGCTTAG- 
GTAACACCAGAGCATTGTCTAACAACATC-3’ for PCR 
amplification. The expression constructs were transformed 
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into E. coli BL21(DE3) (Invitrogen). A single colony of the 
transformant was grown at 37 °C to an ODeoo of 0.7 and 
then transferred to 23 °C. The cells continued to grow to 
ODeo0 = 1.0. The expression was induced for 6 h with 0.4 
mM IPTG. The cells were harvested by centrifugation using 
JLA 8.1000 rotor in Beckman Avanti J-20 XP centrifuge at 
9400g and stored at —80 °C. 

Purification of 3CL?’’ A cell pellet was thawed briefly at 
room temperature and resuspended in 25 mM Tris (pH 8.0), 
50 mM NaCl, 0.1 mM EDTA, and 5 mM DTT. Subsequent 
steps were carried at 4 °C. Cells were lysed by microflu- 
idization and cleared by ultracentrifugation (Beckman L8- 
M) at 125171 x g for 45 min The cleared lysate was loaded 
onto a 10 mL Q-sepharose high-performance column (Am- 
ersham, Piscataway, NJ). 3CL?" was collected in the flow- 
through and loaded to a 10 mL phenyl sepharose high- 
performance column (Amersham) after adjusting the salt 
concentration to 0.6 M (NH4)2SO, and 0.5 M NaCl and pH 
to 7.5. The column was washed with 15—30 column volumes 
of 25 mM Tris-HCl, pH 7.5, 0.5 M NaCl, 0.6 M (NH4)2- 
SO,, 0.1 mM EDTA, and 5 mM DTT. 3CL?° was eluted 
with a 15-column-volume linear gradient to 25 mM Tris, 
pH 7.5, and 5 mM DTT. Fractions containing 3CL"° were 
identified by SDS—PAGE, pooled, and dialyzed overnight 
against 25 mM Tris-HCl, pH 7.5, 100 mM NaCl, 0.1 mM 
EDTA, and 5 mM DTT. Concentrated 3CL? was loaded 
onto a HiPrep 26/60 Sephacryl 200 column (Amersham) 
equilibrated in the dialysis buffer, and eluted in the same 
buffer. 3CL? was pooled, concentrated to 10-20 mg/mL, 
flash frozen, and stored at —80 °C. Protein concentration 
was measured spectrophotometrically using calculated molar 
extinction coefficient of 32590 M7! cm7!. Micromass Q-TOF 
microanalysis gave a single mass peak at 33844 Da, in 
agreement with the calculated mass of 33845 Da for the 
mature 3CL””, indicating that the initiation methionine was 
completely removed. 

Isothermal Titration Calorimetry (ITC). Isothermal titration 
calorimetry experiments were performed using a VP-ITC 
titration microcalorimeter from Microcal LLC (Northampton, 
MA). The syringe was filled with a solution containing 186.5 
UM 3CL?", and the cell was filled with a solution containing 
10 uM inhibitor. Both samples were in 20 mM HEPES (pH 
7.0), 1.0 mM EDTA, 0.005% Tween 20 with 4% DMSO 
added and were degassed immediately prior to use. The 
protein had been exhaustively dialyzed against 20 mM 
HEPES (pH 7.0), 1.0 mM EDTA, 0.005% Tween 20 at 
4 °C. The inhibitor solution was prepared by diluting a 10 
mM stock solution in 100% DMSO with the protein 
dialysate. The final DMSO concentration was 4%. To 
minimize the heats of dilution, 4% DMSO was added to the 
dialyzed protein solution. The titrations were performed at 
25 °C. Ten minutes after a preliminary 2 wL injection a series 
of twelve 10 uwL injections spaced 60 min apart were 
performed. In order to make appropriate corrections for heats 
of dilution, a machine blank (buffer into buffer), a macro- 
molecule dilution control (protein into buffer), and a ligand 
dilution control (buffer into inhibitor) were performed. After 
correction for the heats of dilution, the data were fit using 
the ORIGIN software package provided with the instrument. 
The data were fit using a simple one-set-of-sites model. 
Because the binding was irreversible, only the fitted stoi- 
chiometry was reported. 
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3CL?”? HPLC Assay. The 3CL?"® HPLC assay measured 
the protease-catalyzed cleavage of unlabeled peptides (pep- 
tides 1 and 2). Proteolytic activity of 3CL° was measured 
over 5 min at 25 °C in a 50 uL reaction containing 200 nM 
3CL?®, 20 mM HEPES (pH 7.5), 1 mM EDTA, 200 uM 
substrate, 1.2% DMSO, 0.005% Tween-20, and 2 mM DTT. 
The reaction was terminated with the addition of 200 uL of 
acetonitrile/TFA/water (10%/0.1%/89.9%). The peptide ab- 
sorbances (215 nm) were analyzed by reverse-phase HPLC 
(HP 1100 HPLC system) on a Zorbax SB-C8 column (2.1 
mm x 50 mm, Agilent). Peptides were resolved using a 20 
min 10—95% linear gradient of acetonitrile in 0.1% TFA. 
For example, peptide 1 (SITSAVLQSGFRKMA) was cleaved 
to SITSAVLQ and SGFRKMA. Both products were isolated 
and characterized by mass spectrometry to confirm their 
identity. Chemically synthesized peptide products were used 
as standards. 

3CL?” FRET Assay. For enhanced sensitivity and continu- 
ous monitoring of the proteolytic activity of 3CL?®, continu- 
ous fluorescence resonance energy transfer (FRET) assays 
were developed based on peptides 3 and 4. The 3CL° FRET 
assay measures the protease-catalyzed cleavage of peptide 
3 Tamra-SITSA VLQSGFRKMAK-(Dabcyl)-OH to Tamra- 
SITSAVLQ and SGFRKMAK-(Dabcyl)-OH. Intial rates 
were measured by following the fluorescence of the cleaved 
Tamra (ex 558 nm, em 581 nm) peptide using a Tecan Safire 
fluorescence plate reader at room temperature over the course 
of 10 min. Typical reaction solutions contained 20 mM 
HEPES (pH 7.0), 1 mM EDTA, 4.0 “uM FRET substrate, 
4% DMSO, 0.005% Tween-20, and 2 mM DTT. A typical 
signal was 32000 RFU with a background of 7400 RFU. 
Assays were initiated with the addition of 25 nM 3CL?®. 
The extent of the reaction was determined by producing a 
standard curve with chemically synthesized peptide fragments 
(SynPep Dublin, CA). The experiments were performed 
under first-order kinetics by maintaining 3CL"° substrate 
turnover at < 10%. A second FRET substrate, peptide 4 was 
used to characterize proteolytic activity of 3CLP'® (ex 355 
nm, em 558 nm). Typical reaction solutions were the same 
as above with the Tamra—Dabcy] substrate with the excep- 
tion of the buffer where 20 mM BTP (pH 7.0) was substituted 
for HEPES. A standard curve with the pair of chemically 
synthesized Tamra- and Dabcyl-labeled cleavage products 
(Pfizer, La Jolla) was generated to determine the amount of 
SARS-3CL?°-mediated turnover in a typical FRET assay. 
The amount of cleavage was determined to be under 10%. 
We evaluated the enzymatic activity as a function of enzyme 
concentration with 20 uM peptide 3 and found a linear 
relationship (R = 0.998) between 10 and 300 nM, indicating 
that the assay accurately measured product formation. 
Pseudo-first-order conditions (4 “M peptide 3, 200 nM 
3CL"°) were used to assess the substrate specificity (Keat/ 
K,,) for peptide 3. Under pseudo-first-order conditions, the 
first-order rate constant (k) that describes the full progress 
curve when divided by the enzyme concentration will yield 
the kea/Km value. For peptide 3, the first-order rate constant 
k = 0.00280 + 0.00002 s~! obtained with 200 nM 3CL?° 
yielded a kea/Km value of 14000 + 88 M7! s7!. 

3CL?”’ Colorimetric Assays. The colorimetric assay amide 
substrate, peptide 5, SITSAVLQ-pNA, was used to further 
characterize the proteolytic activity of 3CL?. A typical 
reaction of 0.2 mL contained 12 mM BTP (pH 7.0), 4% 
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DMSO, 2 mM DTT, and 0.005% Tween-20. The absorbance 
produced by the cleavage of pNA (400 nm) was measured 
in cuvettes using a Beckman DU650 spectrophotometer at 
room temperature over the course of 10 min. The amount 
of pNA released was calculated using the molar absorptivity. 
This molar absorptivity of pNA was determined to be 13200 
M~! cm! under our experimental conditions using analytical 
grade pNA. To ensure that experiments were performed 
under first-order conditions, 3CL? substrate turnover was 
determined to be <10%. 

The colorimetric assay substrate SITSAVLQ-pNP (peptide 
6) was used to characterize the enzymatic activity of 3CL?° 
with an ester substrate. Typical reaction solutions contained 
20 mM HEPES (pH 7.0), 4% DMSO, and 0.005% Tween- 
20. The absorbance produced by the cleavage of pNP (400 
nm) was measured using a Beckman DU650 spectropho- 
tometer at room temperature over the course of 10 min. The 
amount of pNP released was calculated using the molar 
absorptivity 7160 M7! cm7! determined under our experi- 
mental conditions using spectrophotometric grade pNP. 
3CL?" substrate turnover was determined to be <10% to 
ensure the experiments were performed under first-order 
conditions. 

Pre-Steady-State Analysis. Pre-steady-state experiments 
were performed on an Applied Photophysics SX.18MV-R 
stopped flow spectrophotometer fitted with a 20 uL cell at 
25 °C with a dead time of 1.0 ms. The cleavage of peptides 
was monitored either by fluorescence of Edans peptide (ex 
335 nm, em 530 nm) or by absorbance at 400 nm for pNA 
and pNP. Enzyme and substrates were diluted 2-fold in assay 
buffer. Concentrations indicated in the text of figure legends 
are the final concentrations after dilutions. Data were 
collected from the average of at least four separate measure- 
ments. The time course of fluorescence or absorbance can 
be described by a single-exponential equation with a steady- 
state term (eq 1) describe previously (24): 


Y=Be"+Att+ce (1) 


Y is the signal (fluorescence or absorbance) at time ¢. B is 
the amplitude of the burst, and D is the observed first-order 
rate constant governing the burst. A is the linear rate of 
increase in fluorescence or absorbance, and c is the back- 
ground. The data were analyzed by simulation of the time 
courses using Pro-K global analysis/simulation program 
(Applied Photophysics). 

pH—Rate Studies. Studies were performed on the effect 
of pH on both the stability of the protein and the catalytic 
rate. The pH stability profile of 3CL?° protease was 
performed to ensure that velocity changes were not due to 
pH-dependent irreversible denaturation. A 500 nM stock of 
3CL?® was incubated for 10 min in 10 mM MES/BTP buffer. 
An overlapping buffer scheme at constant ionic strength was 
used to avoid buffer artifacts: MES from 5.5 to 6.5 and BTP 
from 6.5 to 9.5. The FRET assay with peptide 3 was used 
to monitor protease activity from pH 5 to 9 in 0.5 pH units. 
The concentrated solutions were diluted 40-fold into a 20 
mM HEPES (pH 7.0) solution, and the activity was 
measured. In addition, the pH profile of Tamra-peptide 
fluorescence was performed to ensure that fluorescence did 
not change with pH. There was no significant difference in 
fluorescence of the fluorescent peptides over the pH range 
of 5 to 9. 
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The kea/Km was measured from pH 5 to 9 in 0.5 pH 
increments. An overlapping buffer scheme at constant ionic 
strength was used to avoid buffer artifacts: 30 mM MES 
from 5.5 to 6.5 and 30 mM BTP from 6.5 to 9.5. Full reaction 
time courses were taken under pseudo-first-order condi- 
tions: 4.0 uM FRET substrate, which is much less than the 
K, value. The progress curves were fit to eq 2: 


[S] = [S,]e “ (2) 


[S] is the substrate remaining at time ¢, [S]o is the initial 
substrate concentration, and k = (Kea/Km)[E]. The log Keat/ 
Km was plotted vs pH. The data were fit to eq 3 for two 
ionizable residues, where c is the maximum value of the Kea/ 
Ky parameter. 


log(y) = log(c/(1 + [H*V/K, + K,/[H"))) (3) 


Measurement of Solvent Isotope Effects. Full progress 
curves were measured under first-order conditions for 3CL?° 
in the presence of H,O or D,O. A 7.2% sucrose viscosity 
control was used to assess the impact of the more viscous 
D,0 (7! = 1.24). The viscosities were measured in triplicate 
at 21 °C with an Oswald viscometer (25). The appropriate 
pH meter correction was used to prepare the buffers (pD = 
meter reading + 0.4). The FRET assay format was modified 
to study the effects of DxO on 3CL"”. [7H]-HEPES and -BTP 
were prepared by at least three rounds of solvating HEPES 
and BTP in D2O and lyophilizing the solution. 3CL"° 
reaction solutions contained 20 mM HEPES or BTP (pH 7.0), 
FRET substrate, 4% DMSO, 2 mM DTT, and 0.005% 
Tween-20. The amount of pNA released was calculated using 
the molar absorptivity 13200 M7! cm! in H,O and 12300 
M7! cm™! in D.O. The amount of pNP released was 
calculated using the molar absorptivity 7160 M7! cm“! in 
HO and 4160 M7! cm™! in D,O determined under our 
experimental conditions using spectrophotometric grade pNP 
(Aldrich Chemical Co.). The kinetic data were fit to eq 1 to 
obtain first-order rate constants. Division of the rate constant 
by the enzyme concentration resulted in kea/Km values. 

3CL?” Proton Inventory Studies. Proton inventory studies 
were performed using either FRET (peptide 4) or colorimetric 
(peptide 5) assay formats. The buffer [7H]-BTP was prepared 
by solvating BTP in D2O and lyophilizing the solution three 
times. Solutions of different H2O/D2O fractions (n) were 
prepared by mixing appropriate quantities of BTP buffers 
made up in H20 and D,2O (pH 7.0, pD 6.6). Nine different 
values of n ranging from 0 to 1.0 were used for each substrate 
and six substrate concentrations were used per n value to 
measure the kinetic parameters. 

Analysis of Proton Inventory Data. The experimental data 
were fit to eqs 4—6 by nonlinear regression (2/, 26). 


k, = kl — n+ no — n+ no] (4) 

k, = kg —n+n¢") (5) 

k, =k / —n+no% (6) 

The reaction parameters are expressed as k,(n) functions of 


deuterium atom fraction 1 present in the isotopic solvent, 
where ky is the reaction parameter in HO and ¢7 and @% are 
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the isotopic fraction factors of the transition-state proton and 
the ground-state proton, respectively. 

Inactivation Studies. The inactivation of 3CL?° by the 
Michael acceptor irreversible inhibitor PF-596162 was 
measured as function of time and concentration to determine 
the inactivation constants. Assay conditions for the inactiva- 
tion of 3CL?"® with PF-596162 contained the following: 25 
nM 3CL?, 4 uM peptide 3, HEPES (pH 7.0), and 0.005% 
Tween-20. PF-596162 concentration was varied from 40 to 
0 «M with nine 2-fold dilutions. Background rates were 
subtracted from the rate data. The progress curves for 
irreversible inhibitors were fit to eq 7 where A equals the 
fluorescence of fully cleaved peptide substrate. Data was 
analyzed with the nonlinear regression analysis program 
X1fit3.05 (ID-BS, Guildford, UK). 


F(t) = (A) — ee *+) (7) 


The k.,; was the pseudo-first-order rate constant for the 
inactivation reaction. The slope (Kops/[inactivator]) of a plot 
of kobs VS inactivator constant is the pseudo-second-order 
inhibition constant, which measures potency of the inhibitor 
(880 M7! s~!). Alternatively, the kop; was converted to tip 
and plotted as a function of the inverse of inhibitor 
concentration to yield a Kj = 7.8 uM with a Kinace = 0.012 
s_!. Ki describes the binding equilibrium, and Kinact describes 
the chemical rate of inactivation. 

The pH-dependent inactivation of 3CL**° (100 nM) by 
iodoacetamide (IA) was measured in a reaction mixture 
containing 20 mM BTP—MES buffer combination from pH 
5 to 9 in 0.5 pH units. The FRET substrate peptide 4 was 
used at a 25 uM concentration with 0.001% Tween-20 and 
4% DMSO. The iodoacetamide concentration was varied at 
3, 1, and 0.1 mM final concentration. Inactivation progress 
curves were fit to eq 7 to obtain the first-order inactivation 
constants (k.»s) The second-order rate constant (k,»/[[A]) was 
calculated by dividing the first-order rate constant (Kops) by 
the concentration of iodoacetamide. Data were fit with 
Kaleidagraph 3.5 (Synergy Software, Reading, PA) to the 
equations with pH-independent rate constants and pK, values 
of catalytically competent functional groups. For two rate 
constants and two pK, values, eq 8 was used. For one rate 
constant and one pK, value, eq 9 was used. 


Kopf = kL + 10° PH) 4 4 QPHoPK2)y) 4 
k [1 + 10°%2-P)) (8) 


Kopf = k[1/ + 10° Py] (9) 


where k; and kz are pH-independent rate constants and pK; 
and pK» stand for the pK, values of the catalytically 
competent functional groups. 


RESULTS 


Enzyme Assay and Characterization. Since the foundation 
of meaningful mechanistic studies is based on highly purified 
enzyme, we sought to characterize the fraction of 3CL”® that 
was active using an irreversible substrate mimetic (PF- 
596162). PF-596162 is a peptidomimetic Michael acceptor 
containing a weakly vinylogous ester with the following 
composition: 4-(2-[(4-methoxy-1H-indole-2-carbonyl)-amino]- 
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FicurE 1: Structure of 3CL?® inactivator PF-596162. 


4-methyl-pentanoylamino)-5-(2-oxo-pyrrolidin-3-yl)-pent-2- 
enoic acid ethyl ester (Figure 1). 3CL"® was inactivated by 
PF-596162 as described by the following kinetic param- 
eters: Kj = 7.8 UM; Kinact = 0.012 s~! (Figure 2A,B). 
Crystallographic studies have shown that the inhibitor forms 
a single covalent complex with the active site cysteine residue 
Cys145 (Jay Davies, unpublished results). A similar irrevers- 
ible peptide mimetic inactivator has been shown to form a 
single covalent adduct with Cys145 of 3CL?" (27). Isother- 
mal titration calorimetry with PF-596162 was used to support 
the existence of a single binding site per monomer and to 
determine the fraction of the active sites occupied (Figure 
2C,D). The raw titration data obtained using 60 min injection 
intervals was used to ensure that the irreversible covalent 
binding reaction was allowed to go to completion prior to 
making a subsequent injection (Figure 2C). The data were 
integrated after correcting for heats of dilution (Figure 2D). 
A stoichiometry of 0.96 + 0.01 (moles of ligand per mole 
of 3CL? monomer equivalents) was obtained after fitting 
the data to a one-set-of-sites model. Calorimetry studies 
coupled with the pre-steady-state studies (see below) dem- 
onstrated that 3CL?° was highly pure and fully active. 


To measure protease activity for both amide and ester 
substrates, multiple assay formats (HPLC, FRET, and 
colorimetric) were developed. 3CL?"°-mediated cleavage of 
an untagged peptide was monitored by HPLC for the 
generation of two peptide products and the loss of substrate. 
Based on the conservation of coronavirus substrate specifici- 
ties (28) and modeling studies (7), two substrates for the 
HPLC peptide cleavage assays (peptides 1 and 2) were 
designed. 3CL?"° cleaved peptide | (SITSAVLQSGFRKMA) 
more effectively than peptide 2 (SITSAVLQSGFRK): 45% 
vs 28% cleavage. Kinetic studies in the HPLC assay format 
were performed on the more active peptide, peptide 1, and 
demonstrated that it was a good substrate (Kea/Km = 14000 
M7! s~!) (Table 1). Although the HPLC provided a direct 
measurement of protease activity, its throughput and sensi- 
tivity limited its utility for some mechanistic studies. A FRET 
assay was developed based on peptide 3, which was an 
amino-terminally Tamra-tagged, carboxy-terminally Dabcy]- 
tagged version of peptide 1. While a sensitive tool for 
monitoring activity, peptide 3 had low solubility (25 uM, 
pH 7.0), which prevented the evaluation of saturating 
concentrations. Pseudo-first-order conditions (4 uM peptide 
3, 200 nM 3CL?°) were used to assess the substrate 
specificity (kca/Km) for peptide 3. The steady-state kinetic 
parameters are shown in Table 1. Comparison of the Kea/Km 
values for the untagged peptide | and the tagged peptide 3 
demonstrated that there was minimal interference of the 
FRET tags at positions Ps and Px’ on the 3CL?"°-mediated 
hydrolysis reaction. A second FRET substrate (peptide 4) 
was evaluated based on a previously reported 3CL?” substrate 
(8). Peptide 4 was found to be more soluble (> 100 uM at 
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Ficure 2: Evaluation of 3CL?° inactivator PF-596162: (A) inactivation of 25 nM 3CL?'° with 4 uM peptide 4 by titration with PF- 
596162; (B) plot of ti. vs the inverse of the inhibitor concentration used to calculate the Kj and Kinact for PF-596162; (C) isothermal titration 
calorimetry of 186.5 uM 3CL?" titrated into 10.0 uM PF-596162 at 25.0 °C using 60 min injection intervals; (D) binding isotherm corrected 
for heats of dilution and fit to a one-set-of-sites model. 


Table 1: The Steady-State Kinetic Parameters for SARS 3CL?° Peptide Substrates* 


saturation kinetics pseudo-first-order kinetics 


peptide peptidic substrates Kea (S~!) Kn (uM) Keal Km (M7! s~!) Keal Km (M7! s~!) 
1 SITSAVLQSGFRKMA 8.5 + 0.36 600 + 60 14000 + 1800 b 
3 Tamra-SITSAVLQSGFRKMAK-Dabcyl b b b 14000 + 90 
4 Dabcyl-KTSAVLQSGFRKME-Edans 1.5+0.2 4545 34000 + 7000 34000 + 1500 
5 SITSAVLQ-pNA (amide) 0.86 + 0.06 180 + 30 4800 + 200 5000 + 300 
6 SITSAVLQ-pNP (ester) 0.60 + 0.03 7348 8300 + 1200 b 


“ Peptide | was measured in an HPLC format. Peptides 3 and 4 were measured in a FRET format. Peptides 5 and 6 were measured in colorimetric 
formats. Due to poor solubility of peptide 3, kcat and Km could not be determined. There is a good correlation between the kinetic data obtained from 
full substrate profiles compared with that from pseudo-first-order kinetics (full progress curves). ’ Not determined. 


pH 7.0) and an excellent substrate (Table 1). As measured 
by Keat/Km, peptide 4 was 2-fold better than peptide 1 or 3. 
Peptide 4 is very similar to peptide 3 with the first points of 
divergence from the scissile bond outward at P7 and P7’. A 
small increase in the cleavage reaction also was observed 
for peptide 1 compared with peptide 2, which may be due 
to the additional P¢’ and P7’ positions. To evaluate the effect 
of amide hydrolysis with a good leaving group, a colorimetric 


assay was developed to measure the hydrolysis of peptide 5 
to release pNA. This peptide has the identical P,—Ps 
sequence as peptides 1 and 3. With our 3CL?° protein 
construct, peptide 5 was a good substrate but had 3-fold less 
activity as measured by kea/Km compared with peptides | 
and 3 with their natural P’ amino acid sequence. Finally, to 
study ester hydrolysis, an assay to measure the release of 
pNP from peptide 6 was developed. Peptide 6 was designed 
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FIGURE 3: Pre-steady-state analysis of the solvent isotope effect 
on hydrolysis of peptide 4. The 3CL?'° pre-steady-state reaction 
time course in H,O buffer (A) or D.O buffer (O). Enzymatic 
reactions were performed in 12 mM HEPES (pH 7.0/pD 6.6 
mixtures) with 1 uM 3CL?" and 12.5 uM peptide 4. 


to be identical to peptide 5 except that it had an ester as the 
scissile bond. The measured chemical reactivity (100 uM 
peptidic substrate, pH 7.0, 22 °C) of peptide 6 to nonenzy- 
matic hydrolysis was shown to be very large compared with 
peptide 5. The nonenzymatic rate of peptide 6 (5.0 pmol/s) 
was 3 orders of magnitude faster than that for peptide 5 
(0.0026 pmol/s). Interestingly, by comparison of Kea/Km 
values, the ester was only a marginally better substrate 
compared with the highly related peptide 5 (Table 1). To 
confirm the kinetic findings and to provide methods ame- 
nable for additional studies, kea/Km values were determined 
under pseudo-first-order conditions (Table 1). The kea/Km 
values were invariant of the kinetic method used to determine 
them. 

Pre-Steady-State Analysis. Pre-steady-state kinetic studies 
were initiated in order to better understand the acylation and 
deacylation reactions. Initial pre-steady-state studies of 3CL?° 
showed biphasic progress curves, which are indicative of a 
double displacement (ping pong) mechanism (Figure 3). For 
our studies, we began by using a simple double displacement 
kinetic mechanism as a model system (Scheme 1). Pre- 
steady-state studies of protease-mediated ester hydrolysis by 
a simple double displacement mechanism described by 
Scheme | have been previously described (24). Three main 
parameters were obtained by pre-steady-state analysis: the 
burst rate (b), the burst amplitude (B), and the pre-steady- 
state slope (A). In the biphasic profile, b describes the rate 
of the initial phase, B describes the amplitude of the first 
phase, and A describes the rate of the second phase. The 
pre-steady-state slope should correspond to the steady-state 
rate. 

Because of the high solubility of peptide 6 and its relatively 
low Ky, value, substrate saturation of the ester substrate 
(peptide 6) was possible in the pre-steady-state. The pre- 
steady-state behavior of peptide 6 was studied at multiple 
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FIGURE 4: 3CLP° pre-steady-state evaluation of pre-steady-state 
slope and burst amplitude as a function of substrate concentration: 
(A) reciprocal plot of the pre-steady-state slope as a function of 
peptide 6 concentration; (B) reciprocal plot of square root of the 
burst amplitude as a function of peptide 6 concentration. 


concentrations (12.5, 25, 50, 100, and 200 uM). The burst 
rate (b) was independent of the peptide 6 concentration, while 
the burst amplitude (B) and the pre-steady-state slope (A) 
exhibited saturation kinetics. The pre-steady-state data for 
burst amplitude (B) was plotted as a function of peptide 6 
concentration to the following equation: B~°°= ((k2 + k3)/ 
ky)E°? + ([Kimapp(k2 + k3)]/k2)SE°°] (Figure 4B) (24). The 
approximation of k, > k; allowed for the estimation of the 
enzyme concentration from the y-intercept, 2.2 + 0.3 uM, 
which matches the estimated added protein concentration (2.0 
uM). At the maximum concentration of peptide 6 (200 uM), 
the burst amplitude (B) was determined to be 2.0 + 0.2 uM, 
which also supports these approximations. From the deter- 
mined y-intercept and the known enzyme concentration (2 
UM), the value of k/(k2 + k3) was determined to be 1.0 + 
0.1. To achieve a value close to unity, k2 must be much 
greater than k3. For the mechanism described in Scheme 1, 
the burst amplitude (B) simplifies to B = [E]o when kz > ks 
and [S]o >> Kmapp- For 3CL?", the maximum burst amplitude 
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Table 2: Steady-State Analysis of the Solvent Isotope Effects on Substrate Hydrolysis“ 


H,0 D20 solvent isotope effect 
peptide Keat (s—!) Kn (uM) Kea! Km (M7! s) Keat (s“!) Kn (uM) Kea! Km (M7! s!) Kea) Ky) Keav Ky!) 
3 b b 14000 + 90 b b 24000 + 1600 b b 0.58 + 0.040 
4 1540.2 4545 34000 + 7000 3.7:2'0,1 4447 90000 + 11000 0.39 + 0.05 1040.2 0.37+0.06 
5 0.86 +0.05 180+ 30 4800 + 200 0.98+ 0.03 175+ 16 5600 + 200 0.88 + 0.07 1040.2 0.83 + 0.06 
6 0.60+0.03 7348 8300+ 1200 0.47+0.03 7447 6500 + 1200 1.30+0.14 0.99 + 0.33 1.3 + 0.4 


“ Kinetic constants for peptides 4, 5, and 6 were measured from variable substrate concentration studies, while that for peptide 3 was measured 
under pseudo-first-order conditions due to insufficient solubility to achieve saturation. ’ Not determined. 


Table 3: Pre-Steady-State Analysis of 3CL?° Peptide Hydrolysis of Amide and Ester Substrates Including Solvent Isotope Effects“ 


H,0 D20 solvent isotope effect 
peptide B (uM) b(s"!) A (pmol/s) B (uM) b(s"!) A (pmol/s) B (uM) b(s"!) A (pmol/s) 
4 0.82 + 0.02 5042 26 + 0.3 1.1 + 0.06 28 + 0.6 47+ 0.6 0.76 + 0.04 1.8+0.1 0.55 + 0.01 
5 1.4+0.3 23 40.5 1441.3 1.2+0.3 2442 1741 1.1403 10+0.1 0.84 + 0.07 
6 1.6 £ 0.08 23 £0.8 10.9 + 0.3 0.95 + 0.11 17.2 2.2 7340.4 0.61 + 0.07 1.3+40.2 1.5 + 0.03 


« The amplitude of the burst (B), first-order rate of the burst (b), and steady-state rate (A) were measured. Peptide 4 was evaluated at 12.5 uM 
(Km = 44 «M) with 1 “M 3CL?". Peptide 5 was evaluated at 200 uM (Km = 180 uM) with 1 uM 3CL”®. Peptide 6 was evaluated with 100 uM 


(Km = 73 uM) peptide and 2 uM 3CL?"°. 


(B) was equivalent to the amount of enzyme added. If 
Scheme | holds, this indicates that kz > k3 and [S]o => Kmsapp- 
The pre-steady-state slope data were fit to the following 
equation: 1/A = (W/key)JE + Kmyapp/(KcaS) (24) to determine 
Keat = 0.43 + 0.05 s~! and Knapp = 73 = 7 UM (Figure 4A). 
The ka, and Ky» values calculated from the evaluation of the 
pre-steady-state slope (A) were very similar to the values 
determined through steady-state analysis (Table 1). For a 
simple double displacement mechanism (Scheme 1), Kear = 
kok3l(ko + ks). If ko > ka, then Keat simplifies to k3. From the 
maximum pre-steady-state slope of peptide 6, we can 
estimate k3 to be 0.43 s~!. Finally, the burst rate data were 
analyzed. The burst rate b = ka[S]o/(Ks + [S]o) where K; is 
the substrate binding constant. If ky > ks, and [S] > Ks, 
then the burst rate (b) is equal to the acylation rate (kz). With 
K, > [S]o, the equation simplifies to b = ka[S]o/K;. From 
burst rate of peptide 6, we can estimate the lower limit of 
the acylation rate (k2) to be 18 s~!. The relationship Knapp 
= Kk3/(k2 + kz) has been previously described (Scheme 1) 
(24). With the determined value of ko, k3, and Kmapp, We Can 
calculate the K, to be 2900 uM. K, can be defined for Scheme 
1 as (kK_-1 + ko)/k,. If there is a high commitment to catalysis 
(kz > k_1), then K; can be simplified to ko/k;. With estimates 
of K, and ko, the on-rate (k,) can be estimated to be 6200 
M7! s~!. Another way to determine k, is from the kea/Kmn 
value (Kea/Km = kiko/[k-1 + ko]). When there is a high 
commitment to catalysis (kz > k-;) then Kea/Km = ki (29). 
The steady-state kea/Km value is similar to the k; determined 
with pre-steady-state data. This correlation indicates that the 
estimates of k2, k3, and K, are reasonable. The estimate of k; 
is substantially lower than would be expected for simple 
diffusion-controlled binding of a substrate to an enzyme, 
which indicates that the process may be more complex. 
With the knowledge from the evaluation of the ester 
substrate, peptide 6, we tried to gain mechanistic insight by 
studying the related amide substrates (peptides 4 and 5). 
Unlike the ester substrates, the amide substrates have a 
hydrogen bond donor in the scissile bond, which has been 
shown to have a significant role in catalysis for other 
proteases (23, 30—33). Because of insufficient solubility and 
higher K,, values, peptides 4 and 5 could not be saturated 


under pre-steady-state experimental conditions. Nevertheless, 
mechanistic insight can be gleaned from the evaluation of 
the amide substrates (peptides 4 and 5) relative to the ester 
substrate (peptide 6). The three substrates were tested at Ky, 
levels. If kz > k3, and [S] => Kmapp, then A = kea[E]o. The 
first condition is met (kz > k3) because biphasic time-course 
profiles were observed for all three peptides, but the second 
condition was not met. The steady-state k.a values can be 
compared with the pre-steady-state velocities by the follow- 
ing relationship: A/[E]o = Kea (Tables 1 and 2). The 
correlations between the steady-state ka values and pre- 
steady-state slope values were good even though the pre- 
steady-state substrate concentrations were at Kn levels. 
Because the velocity at Km should be only 2-fold different 
than that at Vinax, the similarity of the k.4 values derived from 
steady-state and pre-steady-state data is not unexpected. The 
amplitudes of the bursts (B) were similar to the concentration 
of added enzyme for peptides 4 and 5 (Table 3), even at 
apparently subsaturating substrate levels. The burst rate (b) 
of peptide 4 was 2-fold larger than that of peptide 5 or 6 
(Table 3). The burst rates of peptides 5 and 6 were equivalent 
even though peptide 5 has an amide as the scissile bond while 
peptide 6 has an ester bond. Again, the enzyme-catalyzed 
hydrolysis of the more chemically labile ester substrate was 
not substantially more efficient than the corresponding amide 
substrates. The steady-state and pre-steady-state analysis of 
ester hydrolysis provided estimates of the rate constants for 
the individual steps of the kinetic mechanism. Pre-steady- 
state analysis of the amide substrates raised doubts that 
3CL?" has a simple double displacement kinetic mechanism 
because of the similarity of the ester and amide burst rates 
and pre-steady-state slopes. 

Solvent Isotope Effect. Solvent isotope effect (SIE) studies 
are mechanistic tools useful in characterizing catalytically 
relevant proton-transfer reactions, substrate interactions, and 
conformational changes. SIEs probe the impact of D,O on 
the enzymatic kinetic parameters. But since D2O may have 
many effects (enzyme stability, conformational changes, 
substrate binding, and proton transfers in catalytic steps), 
care needs to be exercised in the interpretation of SIEs (26). 
For example, the higher viscosity of D.O could affect a 
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physical step of the reaction. A sucrose microviscogen with 
the viscosity of D2O (7! = 1.24) was determined to have 
no significant effect on the steady-state 3CL? reaction. 
Analyses of the SIEs for peptides 4, 5, and 6 were conducted 
in both steady-state and pre-steady-state. Normal and inverse 
solvent isotope effects were observed in the steady state with 
3CL?° (Table 2). To accurately characterize the SIEs, the 
SIE was determined by both pseudo-first-order conditions 
and a full substrate titration where possible. When pseudo- 
first-order conditions were used to monitor the full turnover 
of 1 uM peptide 4, a significant inverse SIE was observed: 
Keal/Km(H20) = 34000 + 1500 M7! s7!; kea/Km(D20) = 
60000 + 1120 M7! s7!; (Kea/Kin)2029 = 0.57 + 0.04. Under 
similar conditions, the related peptide 3 also displayed a 
significant inverse SIE: peptide 3 (Kea/Km)"20™° = 0.58 + 
0.04. The enhanced solubility of peptide 4 allowed for 
saturation kinetics. With data from a full substrate-concentra- 
tion profile, a large inverse SIE was observed for hydrolysis 
of peptide 4: (Kea/Km)"292° = 0.37 + 0.06. Peptide 5, which 
contained an activated amide as the scissile bond, yielded a 
smaller inverse SIE observed under pseudo-first-order condi- 
tions: Kead/Km(H2O) = 5000 + 300 M7! 87}; Keat/Km(D20) = 
6800 + 200 M7! s7!5 Kea/Kmt202°= 0.73 + 0.06. A full 
substrate titration confirmed the inverse SIE: Keat/Km'#202° 
= 0.83 + 0.06 (Table 2). We determined the SIE for the 
ester substrate, peptide 6, which only differs from peptide 5 
by having oxygen in the scissile bond. In contrast to peptide 
5, there was a normal SIE observed on k. for peptide 6 
(Table 2). Due to the large nonenzymatic rate for peptide 6, 
the analysis of the full progress curve was not attempted. 
Nonetheless, significant differences in the observed SIE 
values for ester and amide substrates were observed. 

Because the 3CL?'° mechanism has both acylation and 
deacylation half-reactions, we sought to study the SIE in the 
pre-steady-state to characterize the acylation reaction (Table 
3). Due to solubility issues encountered with pre-steady-state 
experimental conditions, we investigated the SIE under K,, 
conditions. The differential levels of saturation should affect 
the magnitudes of the effects in the pre-steady-state but not 
the character of the SIE (normal vs inverse). There was a 
large normal SIE on the burst rate of peptide 4 and an inverse 
SIE observed on the pre-steady-state slope (Figure 3). Peptide 
5 did not have an SIE on the burst rate and had a small 
inverse SIE on the pre-steady-state slope. The ester substrate, 
peptide 6, had normal SIE values for both burst rate and 
pre-steady-state slope. The proton-transfer reaction of the 
cysteine sulfhydryl group is unique among amino acid side 
chains because of its low bending and stretching frequencies 
resulting in an inverse SIE (26, 34). The normal SIE on the 
burst rate is not consistent with a rate-limiting general base 
abstraction of a proton from the sulfhydryl group of a 
catalytic cysteine residue. The SIE analysis in the steady- 
state and pre-steady-state indicates that the 3CL? mechanism 
may be more complicated than the minimal protease mech- 
anism (Scheme 1). 

3CL?” Proton Inventory. To gain further understanding 
of the underpinnings of the observed SIEs, the proton 
inventory technique was used to further characterize the 
hydrogenic sites that produced the SIEs. Proton inventories 
study the effect of mixtures of H2O and DO on kinetic 
constants in the steady state (26, 34). The individual 
fractionation factors ¢' and ¢° measure the deuterium 
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FiGurRE 5: The nonlinear, least-squares fit of the kca/Km proton 
inventory data for peptides 4 and 5 (amides). Enzymatic reactions 
were performed in 20 mM HEPES (pH 7.0/pD 6.6 mixtures) with 
400 nM 3CL?" and subsaturating substrate concentrations (5 “M 
peptide 4 or 100 uM peptide 5): (A) data for peptide 4 was fit to 
eq 4, the model for two fractional factors in both the ground state 
and transition state; (B) data for peptide 5 with equivalent fits to 
the equation for fractionation factors in both the transition state 
and ground state (eq 4) or just the ground state (eq 6). 


preference of a particular site (in the ground or transition 
state) relative to the deuterium preference for the bulk water 
(34). Bonds to neutral oxygen or nitrogen are essentially 
equivalent to those in bulk water and often lead to a 
fractionation factor of unity for #°%. The reciprocals of the 
fractionation factors represent the contribution of the site to 
the overall solvent isotope effect. Since the observed SIEs 
for peptides 4 and 5 had little contribution from Ky, we 
performed a proton inventory by monitoring the full turnover 
of a substrate under subsaturating conditions to yield Kea/ 
Kn values. Data for peptide 4 were fit to multiple forms of 
the Gross—Butler equations, eqs 4—6 (26, 34) (Figure 5A). 
The fit of peptide 4 data to the model with fractional factors 
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FIGURE 6: The nonlinear least-squares fit of the k.,, proton inventory 
data for peptide 6 (ester). Enzymatic reactions were performed in 
10 mM HEPES (pH 7.0/pD 6.6 mixtures) with 500 nM 3CL?° and 
a saturating substrate concentration (200 “M peptide 6). The normal 
solvent isotope effect fit best to the equation for one fractionation 
factor in the transition state (eq 5). 


in both the ground state and transition state was better (eq 
4, R? = 0.99) compared with the models for a single fraction 
factor in either the transition state (eq 5, R* = 0.86) or the 
ground state (eq 6, R* = 0.95). The fractionation factors for 
peptide 4 were determined to be $7 = 0.58 + 0.03 and ¢% 
= 0.33 + 0.02. The proton inventory data for peptide 5 were 
fit to eqs 4—6 with equivalent fits to the equation for 
fractionation factors in both the transition state and ground 
state (eq 4, R* = 0.98) or just the ground state (eq 6, R? = 
0.98) (Figure 5B). The fractionation factor for both transition 
state and ground state fractionation factors yielded 6" = 0.95 
+ 0.14 and ¢& = 0.70 + 0.11, while the model for only a 
ground state fractionation factor yielded ¢° = 0.74 + 0.01. 
By either fit, the contribution of @" to the SIE was small. 
With the ester substrate, peptide 6, the observed normal SIE 
resided in ka. Because of the significant background 
hydrolytic rate, the proton inventory was performed under 
pseudo-Vmax conditions (500 nM 3CL?", 200 uM peptide 6, 
3 min reaction times) to obtain an estimate of Kea (Keat**). 
This approach was facilitated by the invariance of the Km as 
a function of D,O content. The proton inventory for peptide 
6 was fit to the equations for either one fractionation factor 
in the transition state (eq 5, 6’ = 0.33 + 0.02, R? = 0.99) 
or fractionation factors in the transition state and ground state 
(eq 4, ¢' = 0.40 + 0.08 and ¢° = 1.1 + 0.2, R? = 0.99) 
(Figure 6). Since the ¢% is unity in the two site model, the 
data supports a @" single site model. 


pH—Rate Profile. In probing the 3CL?° mechanism, we 
sought to define the pK, values of the catalytically relevant 
residues to see whether different substrates required different 
ionizable residues to affect catalysis. With a different 3CL?° 
construct and a peptide (TSAVLQ-pNA) related to peptide 
5, the pK, values have been previously reported to be the 
following: pKj; = 6.25 + 0.04 and pKa2 = 8.29 + 0.04 
(10). Because the current studies were performed with a 
different, more active 3CL? construct, we sought to 
reinvestigate the pH—rate profile and expand the analysis 
to include ester hydrolysis. To ensure valid pH—rate results, 
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FiGurE 7: Alkylation of the catalytic cysteine with iodoacetamide. 
Enzymatic reactions were performed in 10 mM BTP—MES buffers 
(pH 5.5 to 9) with 200 nM 3CL?"° and 50 “uM peptide 4. Different 
iodoacetamide concentrations were used to measure the inactivation 
rates: 3 mM for pH 5.5 to 6.5, 1 mM for pH 7 to 8, and 0.1 mM 
for pH 8.5—9. Inactivation rates were corrected with the no 
iodoacetamide controls. The data fits better to an equation describing 
a complex ionization curve with both bell-shaped and sigmoidal 
terms (eq 8) (—) compared with the fit to a simple dissociation 
equation (eq 9) (***). 


a pH-stability profile was performed on the 3CL?"° enzyme. 
3CL?" was found to be stable to irreversible pH-dependent 
denaturation in the examined pH range. The pH-—rate 
behavior was then characterized for peptide 3 under pseudo- 
first-order conditions to determine the k,)/K,. This method 
of determining Kea/Km was necessary due to solubility 
limitations. Fitting the data for peptide 3 to eq 3 yielded a 
bell-shaped profile with two pK, values (6.2 + 0.1 and 7.7 
+ 0.1). The 3CL?"-mediated hydrolysis of an amide substrate 
displayed a bell-shaped profile consistent with multiple 
ionizable residues contributing to catalysis. 


pH-Dependent Alkylation of 3CL?” with Iodoacetamide. 
The profile of cysteine alkylation as a function of pH is a 
method to provide evidence of a thiolate—imidazolium ion 
pair mechanism. A thiolate—imidazolium ion pair would 
allow for enhanced reactivity at pH values substantially less 
than the pK, of a typical cysteine residue. This method has 
been used for other cysteine proteases (/8). Cysteine 
alkylation of 3CL?'® was feasible because the catalytic 
cysteine residue is preferentially nucleophilic. Inactivation 
studies with the cysteine-reactive, weakly electrophilic 
inactivator PF-596162 showed a 1:1 inactivation stoichiom- 
etry. Another irreversible inhibitor has been shown to 
specifically modify the catalytic cysteine residue of 3CL”° 
(27). The pH dependence of the pseudo-second-order rate 
constant of inactivation of 3CL?° (kg,,/[IA]) with iodoac- 
etamide is shown in Figure 7. Inspection of the pH 
inactivation rate profile shows that it does not follow a simple 
dissociation curve (Figure 7). The rate of inactivation of 
3CL?° by iodoacetamide was significant at neutral to slightly 
acidic pH. The data fit better to an equation describing a 
complex ionization curve with both bell-shaped and sigmoi- 
dal terms (eq 8, R? = 0.998) compared with the simple 
dissociation curve (eq 9, R? = 0.994) (Figure 7). The fit to 
eq 8 resulted in the following values: kj = 5.9 + 0.1 M! 
s—!; ky = 0.29 £ 0.18 M7! s7!; pK, = 5.7 + 0.09; pKi2 = 
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7.9 + 0.1. The fit to eq 9 resulted in the following values: 
k=5.7+0.2 M's"; pK, = 7.8 + 0.1. The data fits were 
evaluated by evaluation of the residuals (35). The residual 
sum of the squares for eq 8 was 0.058, while it was 0.16 for 
eq 9. Evaluation of the residuals to the data fits were 
performed by leaving one data point out and using the 
equations to predict the missing data. The simple model (eq 
9) was less predictive at the lower pH range (pH < 7). The 
pH-inactivation profile (6.0 to 8.5 pH range) of PF-596162 
was examined in a similar manner to iodoacetamide. PF- 
596162 was potent in the examined pH range with the Kops/ 
[PF-596162] only reduced by 50% at pH 6.0. Although the 
reactivity to alkylation agents at neutral pH condition is 
suggestive of an active site ion-pair mechanism, it can only 
be considered as corroborative evidence. 


DISCUSSION 


Steady-State Characterization of the 3CL?"” Activity. The 
observed high enzymatic activity achieved with untagged, 
highly purified 3CL?° was a critical prerequisite to the 
mechanistic analysis described in the current study. Although 
other laboratories have reported a similar high enzymatic 
activity (J, 13), there are 3CL”” studies that use a construct 
that has 100-fold lower activity (6, /0). More recent studies 
have shown that extra non-native N-terminal residues have 
a detrimental effect on enzymatic activity (/3). The catalytic 
activity observed for 3CL?” (Table 1) are similar to or greater 
than the activities reported for other proteases (20, 2/, 36). 
As such, the 3CL?" preparation used in this study appears 
to be appropriate for meaningful mechanistic studies. 

A number of studies have explored the substrate require- 
ments of 3CL?". P4—P2’ residues were reported to be 
indispensable for effective cleavage while P3’—P7’ are not 
essential but have an effect on cleavage efficiency (37). An 
additional study of the 3CL?" substrate specificity of P;—Py 
substrate positions found that the P, glutamine can be 
replaced by a histidine residue, the P2 leucine residue is 
obligate, and the P3/P, positions are more tolerant of variation 
(27). Although peptides 5 and 6 have a P1’ mimetic and no 
P2’, they are good substrates. As such, the P2’ position 
appears to be dispensable for catalysis. Peptide substrates 
with an amide in the scissile position (peptides 1 and 3) are 
cleaved 3-fold more efficiently as measured by Kea/Km 
compared with the peptidic substrates with an activated 
amide (peptide 5) or an ester substrate (peptide 6). This 
indicates that the leaving group may have a role in catalysis. 
Studies of papain by Polgar showed a substantial role of the 
leaving group on the substrate efficiency as measured by 
the kca/Km of a large panel of substrates (33). The role of 
the leaving group on catalysis was fortified by using two- 
protonic state reactivity probes of S1—S2 to evaluate papain 
(32). Typically for serine/cysteine proteases, acylation is 
thought to be rate-limiting for amide hydrolysis while 
deacylation is thought to be rate-limiting for esters (18, 38). 
As such, enzyme-mediated hydrolysis of esters is typically 
faster than that with amide substrates. Analysis of 3CL?"° 
substrates demonstrated that the ester substrate peptide 6 had 
kinetic parameters similar to the highly related aniline 
substrate (peptide 5) despite the large difference in inherent 
reactivity. For papain, there was a large difference in Keat of 
pNA amide substrates (X-Phe-Leu-pNA, 0.35 s~') compared 
with the related pNP ester substrate (X-Phe-Leu-pNP, 7.85 
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s-') (21). For the chymotrypsin hydrolysis of acetyl- 
phenylalanine substrates, the ka of the amide was 1500- 
fold smaller than that for the ester (38). The 3CL?"°-mediated 
hydrolysis of highly related ester and amide substrates may 
be distinct from chymotrypsin and papain. To further 
understand the ester and amide hydrolysis reactions, pre- 
steady-state analysis was pursued. 

Pre-Steady-State Analysis. Pre-steady-state evaluation of 
3CL?° provided significant insight into the kinetic mecha- 
nism of 3CL?°-mediated catalysis of ester and amide 
hydrolysis. The evaluation of ester hydrolysis as a function 
of substrate concentration allowed for the estimation of the 
rate constants of the kinetic mechanism and is consistent with 
a double displacement (ping pong) mechanism. Unlike 
related enzymes (23, 38, 39), the burst rate of 3CL?°- 
mediated ester hydrolysis was shown to be independent of 
substrate concentration. The 3CL?" burst rate observations 
can be accounted for if there was a high commitment to 
catalysis: the acylation rate (ky) is much faster than the 
substrate dissociation rate (k_;). From the burst rate of the 
ester substrate, we can estimate the lower limit of the 
acylation rate (18 s~'), which was 10—100 fold higher than 
that for related enzymes (2/). The pre-steady-state studies 
allowed us to estimate the substrate binding constant (Ks) 
and on-rate (k;). The determined K, and k; values for 3CL?° 
were similar to values reported for related enzymes (2/). 
The ester substrate showed that 3CL?" may catalyze the 
acylation more effectively than papain as determined by ky, 
ky, and k_, values. The detailed kinetic analysis of 3CL?"°- 
mediated ester hydrolysis is consistent with a kinetic mech- 
anism similar to Scheme 1. 

As observed for ester hydrolysis, the pre-steady-state 
progress curves for the amide hydrolysis reactions were 
biphasic. As expected, the 3CL?"°-mediated hydrolysis of 
amide substrates had a tight correlation between the pre- 
steady-state slope values and the steady-state k.a, values. With 
3CL?®, the deacylation values estimated from the pre-steady- 
state slopes are similar to each other: peptide 4, 1.3 s~!; 
peptide 5, 0.71 s~!; and peptide 6, 0.47 s~!. The similarity 
of these values could be because the deacylation reaction 
occurs from a similar acyl-enzyme intermediate. But we 
know from SIE studies that there are significant differences 
in the catalytic events that govern the pre-steady-state rate. 
The pre-steady-state slope may include more than the 
deacylation reaction (P; product release). For the serine 
protease chymotrypsin mediated hydrolysis of pNP esters, 
the burst amplitude (B) equaled the added enzyme concentra- 
tion when kz > k3 and Kmapp < [S]o (24). For 3CL?®, the 
first criterion is satisfied while the second is not. Nonetheless, 
for peptides 4 and 5, the burst amplitude was similar to the 
added enzyme concentration, which again indicates an 
underlying mechanistic departure for the amide hydrolysis 
reaction. 

Evaluation of the burst rate data for amide and ester 
hydrolysis reactions was revealing. The burst rate for peptide 
4 was 2-fold faster than burst rate of the peptides with better 
leaving groups: pNA (peptide 5) or pNP (peptide 6). What 
was unexpected is that the burst rates for both peptides 5 
and 6 tested at K, levels were equivalent (Table 3). Peptides 
5 and 6 form the same acyl intermediate with the scissile 
bond of peptide 5 being an amide while the scissile bond of 
peptide 6 is an ester. As such, the acylation rates are expected 
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to differ. If k, > k3, then the burst rate b = ko[S]o/(K, + 
[S]o). If k2 => kz and [S] > K,, then the burst rate (b) is 
equal to the acylation rate (kz). The detailed study of pre- 
steady-state parameters as a function of ester concentration 
shows that this relationship may not hold and probably will 
not for peptides 4 and 5. With K, > [S]o, the relationship 
simplifies to b = [S]o(k2/K,). If both the amide and ester 
hydrolysis reactions occur with a simple kinetic mechanism 
(Scheme 1), then the similar burst rates should be accounted 
for by differing k2/K, ratios. To evaluate the ratio of k2/K, 
for peptides 5 and 6, one can re-evaluate the steady-state 
data. The steady-state K,, values for peptides 5 and 6 differ 
by less than 3-fold. From the relationship Kmapp = Ksk3/(k2 
+ k3) and ky > k3, we can simplify the relationship to Km,app 
= Kks/ko. Since peptides 5 and 6 form the identical acyl 
intermediate, the deacylation rate (k3) by inference should 
be identical. The similarity of the K,, values should be 
propagated from the ratio of K,/k2 if k3 reflects a deacylation 
rate. If kp for the ester is substantially larger than the k2 for 
the corresponding amide, then the K, must be significantly 
larger (poorer binding). Alternatively, the kinetic mechanism 
of the 3CL°-mediated amide hydrolysis may be more 
complicated. The rate constant k3 for amide hydrolysis may 
not describe only the deacylation reaction but may encompass 
additional mechanistic elements (e.g., Scheme 2). His41 
general-base-catalyzed hydrolysis of the highly reactive ester 
substrate is possible but probably unlikely due to the biphasic 
progress curves observed in the pre-steady state. Because 
the kinetic mechanism describes only the sequence of 
catalytic events and not the details of the catalytic mecha- 
nism, other mechanistic techniques were necessary. 

3CL?P” Solvent Isotope Effects. Solvent isotope effect (SIE) 
studies of enzyme reactions provide valuable mechanistic 
information on proton-transfer reactions and conformation 
changes involved in the transition state. For typical proton- 
transfer reactions, an X—D bond requires more energy to 
break than an X—H bond and a normal SIE is observed 
(v?20/30 < 1). A special case is the proton transfer from a 
sulfhydryl group in which an inverse SIE can occur (P2029 
> 1). An inverse isotope effect can be generated by events 
other than sulfhydryl proton-transfer reactions (e.g., confor- 
mational changes, substrate interactions). The inverse SIE 
ON Kea/Km2°”) previously reported for the 3CL?"°-catalyzed 
hydrolysis of the peptide substrate TSAVLQ-pNA was large 
(~0.4), but the enzyme preparation had only 1% of the 
enzymatic activity used in the present study (/0). Since SIE 
determinations are sensitive to many factors, we characterized 
the SIE with a highly active preparation of 3CL?®. In the 
steady state, a combination of pseudo-first-order and substrate 
saturation kinetic studies were used to demonstrate that there 
was a large inverse SIE on the 3CL?°-mediated hydrolysis 
of the amide bonds and a normal SIE for the ester hydrolysis. 
The results were confirmed through the analysis of the SIE 
on the pre-steady-state slopes. In the pre-steady-state studies, 
both ester and amide substrates exhibited a normal SIE on 
the burst rate. This result is not consistent with a rate-limiting 
general base abstraction of a sulfhydryl proton. One possible 
explanation is that the steady-state 3CL”° SIE amide results 
could be accounted for by an enzyme-facilitated release of 
the first product. Release of the amide P; product could be 
facilitated by an interaction with the NH of the leaving group 
or perhaps by a conformational change. Precedent for this 
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proposed mechanism can be found with other cysteine 
proteases. An interaction of backbone carbonyl of papain 
Asp158 with the NH of the leaving group was identified as 
critical for the observed inverse SIE (2/, 23, 30—33). 
Disruption of the ion pair and formation of a new interaction 
of the imidazolium cation with the substrate was postulated 
to be responsible for the inverse SIE (20, 21). Pre-steady- 
state studies of the cysteine protease actinidin, but not papain, 
showed that a postacylation conformational change was 
required to release the alcohol P; product (23). The key 
residue of actinidin that influences the conformational change 
is Asp161, which is in an analogous position to the papain 
residue (Asp158) that modulates the inverse SIE (23). The 
difference between these residues in the pK, was thought to 
be critical to their role as electrostatic modulators. For the 
serine protease thrombin, the inverse SIE observed is 
associated with a conformational change to release P, (40). 
By either mechanism, to affect P; release, the ester cleavage 
reaction cannot be facilitated by an analogous enzyme— 
substrate interaction. Finally, evaluation of the burst rate data 
provides important insight into the 3CL°"? mechanism. To 
understand the normal SIE observed on the burst rate of 
3CL?", detailed studies of papain serve as a useful model to 
explore. Studies of papain family proteases have demon- 
strated that the architecture of the thiolate—imidazolium ion 
pair occurs at pH conditions that do not support catalysis. 
This led to the concept that an “electrostatic switch” 
modulates the thiolate—imidazolium ion pair by hydrogen- 
bonding interactions of an active site aspartate residue with 
the histidine residue of the ion pair (20, 47—43). Involvement 
of this interaction in the acylation reaction could lead to the 
observed normal SIE in the 3CL”” burst. The solvent isotope 
studies illuminated critical similarities and differences be- 
tween 3CL?°-mediated ester and amide hydrolysis. 

Proton Inventory. The proton inventory technique has been 
used to further characterize the SIE. The proton inventory 
describes the number of exchangeable protons in the enzyme 
that are involved in the catalytic mechanism resulting in the 
observed SIE. This technique can characterize the solvent 
isotope effects at each of the contributing exchangeable sites, 
which are quantified with fractionation factors (@) (26, 34). 
Fractionating sites can appear in either reactant (ground) 
states or transition states. A @ value considerably smaller 
than unity in the transition state may be an indicator of a 
“proton in flight” although there are other interpretations 
(e.g., low-barrier hydrogen bonds). In the reactant state, there 
is simply a fractionating site with a stable potential. Such 
fractionation is expected for low-barrier hydrogen bonds, for 
sulfhydryl groups (where the low force constant for the S-H 
bond favors H over D), and for ion pairs such as thiolate— 
imidazolium, where a broad shallow potential (low force 
constant) also favors H over D. A limitation of the technique 
is that all solutions are not unique. For example, a linear 
proton inventory can be experimentally determined even 
though a complex array of proton-transfer reactions can occur 
(34). With these caveats in mind, we explored the proton 
inventories for both the ester and amide hydrolysis reactions. 
The “bowed down” shape of the proton inventory for the 
hydrolysis of peptide 4 by 3CL*° indicates that there are 
two fractionating sites (exchangeable sites) involved in the 
acylation reaction: one in the transition state and one in the 
ground state. The ground-state exchangeable hydrogenic site 
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could be derived from the tautomerization of the reactive 
thiolate—imidazolium ion pair with the neutral thiol— 
imidazole form. The transition-state fractionation factor could 
be attributed to the disruption of the ion pair and the 
formation of a new interaction of the imidazolium cation 
and the substrate. A PI similar to peptide 4 was observed 
for papain-mediated hydrolysis of pNA peptide substrates 
(21). Further evidence of this mechanism for 3CL?° can be 
garnered from the comparison of PIs with amide and 
activated amide substrates. The transition state has a partial 
bond between the catalytic cysteine residue and the substrate 
carbonyl carbon with the double bond character of the C= 
O reduced (2/). The electronic nature of the carbonyl in the 
transition state will be affected by the leaving group. A 
normal amide will require more assistance than an activated 
amide (pNA) because the pNA can stabilize charge through 
resonance. The proton inventory data for peptide 5 was more 
consistent with a fractionating site in the ground state, while 
peptide 4 had fractionating sites in both the transition state 
and ground state. The 3CL”° PI results for peptides 4 and 5 
are consistent with an ion-pair mechanism for the catalysis 
of the acylation reaction. 

From the kinetic evaluation of ester hydrolysis, we have 
shown that k.a conditions should measure the deacylation 
half-reaction. The observed linear PI for 3CL?"® hydrolysis 
of peptide 6 is consistent with a single fractionating site. A 
possible interpretation is that His41 of the 3CL"° catalytic 
dyad is acting as a general base to facilitate activation of 
water in the deacylation second half-reaction as was postu- 
lated for papain (/9). Taken as a whole, the proton inventory 
results for the three peptides studied indicate a mechanism 
with the following components: (1) initiation of an ion-pair 
mechanism in the first half-reaction, (2) enzyme-facilitated 
release of the first product, and (3) general base-catalyzed 
deacylation reaction. 

pH-Dependent Alkylation of 3CLP’ by TIodoacetamide. 
Inactivation studies with PF-596162 and studies in the 
literature (27) show that the 3CL?” catalytic cysteine residue 
was preferentially nucleophilic. Inspection of the pH inac- 
tivation rate profile of iodoacetamide shows that it does not 
follow a simple dissociation curve (Figure 7). There is a 
wealth of studies of cysteine proteases that have similar pH- 
inactivation profiles (18). The iodoacetamide reactivity to 
3CL?" was significant at neutral to slightly acidic pH. A 
cysteine—histidine ion pair would allow for enhanced 
reactivity at pH values substantially less than the pK, of an 
ordinary cysteine. Interestingly, the pK, values of the 
ionizable residues necessary for iodoacetamide inactivation 
are within error of the pK, values of catalytically relevant 
ionizable residues as determined from the pH—rate profiles. 
Although the reactivity to alkylation agents at neutral pH 
condition is suggestive of an active site ion-pair mechanism, 
it can only be considered corroborative evidence. 

A Model for the Mechanism of the 3CL?"’ Emerging from 
the described studies are the details of the 3CL”° kinetic 
steps, which can be assembled into a coherent mechanism. 
A key finding is that the similar 3CL?"® burst rates for ester 
and amide reactions may be attributed to a common step in 
the hydrolysis reactions. We propose that 3CL?" may have 
an “electrostatic trigger” that is rate-limiting and renders the 
burst rate independent of substrate reactivity. Another finding 
is that the 3CL"°-mediated ester and amide hydrolysis show 
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similar catalytic efficiencies while the SIE and PI studies 
reveal mechanistic departures. If there is an enzyme- 
facilitated P; product release due to a critical NH interaction 
as observed for both papain and actinidin, then only amide 
substrates may be able to access this mechanism. The inverse 
SIE and nonlinear PI results for the amide substrates are 
consistent with an ion-pair mechanism. Enhanced reactivity 
toward iodoacetamide at lower pH conditions is also indica- 
tive of an ion-pair mechanism. The normal SIE and linear 
proton inventory for ester hydrolysis is similar to papain and 
may indicate that the deacylation reaction is general base 
catalyzed. The proposed catalytic model for 3CL?°-mediated 
amide hydrolysis is similar to papain, but the studied 
substrates indicate that it may affect catalysis differently. 
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